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The D-brane world is an idea that we are living on the D-brane imbedded in a 10- or 11-dimensional 
spacetime of string theories, aiming at the construction of realistic models from string theories. We 
investigate the cosmological aspects of the D-brane world, focusing on homogeneous anisotropic 
cosmology driven by the dilaton and the NS-NS 2-form field which becomes massive in the presence 
of the D-brane. The dilaton possesses the potential due to the presence of the D-brane, various form 
field fluxes, and the curvature of extra dimensions. In the absence of stabilizing potential, we found 
the attractor solutions for this system which show the overall features of general solutions. In the 
presence of the non- vanishing NS-NS 2-form field, the homogeneous universe expands anisotropically 
while the D-brane term dominates. The isotropy is recovered as the dilaton rolls down and the 
curvature term dominates. With the stabilizing potential for the dilaton, the anisotropy developed 
by the initial NS-NS 2-form field flux is erased as the NS-NS 2-form field begins to oscillate around 
the minimum, forming the B-matter, and the isotropic matter-dominated universe is obtained. 
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I. INTRODUCTION 

The idea of the brane world first appeared as an al- 
ternative solution to the gauge hierarchy problem by al- 
lowing large or warped extra dimensions while we are 
leaving on a 3-brane to avoid various known experimen- 
tal bounds p|. String theories have such structures as 
D-branes on which open strings can reside and thus pro- 
vide a framework where realistic brane world models 
can be constructed. The D-brane world is an idea that 
we are living on the D-brane imbedded in a 10- or 11- 
dimensional spacetime of string theories, aiming at the 
construction of realistic models from string theory. One 
way to see how realistic the D-brane world models can 
be is to look into the cosmological aspects of them. 

Low energy effective theories derived from the NS-NS 
sector of string theories contain the gravity, g^, the dila- 
tion, $, and the NS-NS 2-form field, B^ v . The existence 
of the last degree of freedom leads to intriguing implica- 
tions in string cosmology 0- In four spacetime dimen- 
sions, the massless antisymmetric 2nd-rank tensor field 
is dual to the a pseudo-scalar (axion) field @], and the 
axion-dilaton system is known to develop an unobserved 
anisotropy in our Universe, which can be diluted away at 
late times only in a contracting universe Q. 
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Such a disastrous cosmological situation can be re- 
solved also in string theory which contains another dy- 
namical object called D-brane The gauge invari- 
ance on the D-brane is maintained through the coupling 
of the gauge field strength to the NS-NS 2-form field 
[M IE Of- Then, the effective action derived on the D- 
brane describes the NS-NS 2-form field as a massive and 
self-interacting antisymmetric tensor field. Cosmological 
evolution of such a tensor field has been investigated in 
Ref. p| assuming that the dilaton is fixed to a reasonable 
value. Although the time-dependent magnetic B field 
existing in the early universe develops an anisotropy in 
the universe, it was realized that the matter-like behav- 
ior of the B field (B-matter) ensures a dilution of the 
anisotropy at late times and thus the isotropy is recov- 
ered in reasonable cosmological scenarios 8J. In such 
sense the effect of antisymmetric tensor field on the D- 
brane is distinguished from that of field strength of the 
U(l) gauge field @. 

In this paper, we investigate the cosmological evolution 
of the dilaton-NS-NS 2-form field system in our Universe 
which is assumed to be imbedded in the D-brane. In 
general, the D-brane and the NS-NS 2-form field couple 
to the R-R form fields. Here, we assume a trivial R-R 
background and omit them in our analysis. The usual 
string cosmology with the dilaton suffers from the noto- 
rious runaway problem, which is also troublesome in the 
D-brane universe. In our study, the dilaton obtains two 
exponential potential terms due to the curvature of extra 
dimensions A, and the D-brane tension (the mass term of 
the B-matter) tub- It is interesting to observe that the 
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dilaton can be stabilized for negative A [lj} which, how- 
ever, leads to a contracting universe due to the effective 
negative cosmological constant in our Universe. When A 
is positive, the B-matter dominance will be overturned by 
the A dominance as the dilaton runs away to the negative 
infinity. As a consequence of this, the initial anisotropy 
driven by the B-matter can also be diluted away at late 
times. 

For a realistic low energy effective theory, string theory 
must be endowed with a certain mechanism generating 
an appropriate vacuum expectation value for the dilaton. 
In such a situation, the dilaton is expected to be stabi- 
lized at some stage of the cosmological evolution affecting 
the dynamics of the NS-NS 2-form field. Taking an ex- 
ample of the dilaton stabilization, we will also examine 
the cosmological evolution of the dilaton-NS-NS 2-form 
field system in which the essential features of Ref. 8] are 
reproduced. 

This paper is organized as follows. In Section 2, we 
describe the low energy effective action of the D-brane 
world and the corresponding field equations. In Section 
3, we find analytic and numerical cosmological solutions 
of the dilaton-NS-NS 2-form field system to observe an 
intriguing interplay of the curvature A and the D-brane 
tension tub- In Section 4, we consider the evolution of the 
anisotropic universe with the dilaton stabilization which 
leads to a satisfactory cosmology of the D-brane universe. 
We conclude in Section 5. 



II. EFFECTIVE FIELD THEORY OF THE 
D-BRANE WORLD 

The main idea of the D-brane world is that we reside 
on a Dp-brane imbedded in 10- or 1 1-dimensional space- 
time with extra-dimensions compactified. The bosonic 
NS-NS sector of the D-brane world consists of the U(l) 
gauge field living on the Dp-brane and the bulk de- 
grees including the graviton g^, the dilaton <i>, and the 
antisymmetric tensor field of rank-two B^ v . 

The low energy effective action consists of the bulk 
action 



S_ 



NS 



i? + 4(<9$) 2 - -^H 2 



and the brane action in p spatial dimensions; 
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Recall that a' defines the string scale; m s = a' -1 / 2 . Here 
the spacetime indices are not explicitly expressed and 
the field strength H of the NS-NS form field is H^ p = 
0[ M S„ p ] because of Bianchi identity of the U(l) gauge 



field. In the presence of the brane, the gauge invariance 
of is restored through its coupling to a U(l) gauge 
field and the gauge invariant field strength is [§j 



B u 



B„ 



2na'F fll/ , where = d^A v - d v A p . (3) 



The D-brane and the NS-NS form field on it couple 
to the R-R form field as shown in the second term of 
Eq. J5J. In this paper, we consider the case that the D- 
brane and the NS-NS form field are homogeneous in our 
4D world and all R-R form fields live in extra dimensions 
so that their sole effect is the curvature term in the 4D 
effective action which will be discussed below. Thus all 
R-R form fields are not considered from now on, assuming 
the trivial R-R background. We leave the more general 
case of both the NS-NS and the R-R form field living in 
our 4D world as our future work. 

After the compactification of extra dimensions, the 
four-dimensional effective action of the bosonic sector in 
the string frame is || 



s s = Tr? I d^x^J-g 



-2* 



( /? - 2A + 4V P $V^$ - ^H^ p H^p 



,2 -* 



1 + 



,(4) 



where the tilde denotes the string frame quantity and 
\\J— ge-fiuafsB a ^ with £0123 = 1. The parameter 



B* 



tub is defined by m 2 B = 2k, 2 Hp- If we assume that the 
six extra-dimensions are compactified with a common ra- 

1 o p — 3 15 — p 

dius i? c , one finds m B — ^ (g 2 /4ir)^5- (R C M P ) 8 M P 
where g s = e* and M P = 2.4 x 10 18 GeV is the four- 
dimensional Planck mass. The qualitative features of our 
results do not depend on specific values of p. Thus, we 
take p = 3 for simplicity. The A term comes from the 
scalar curvature or the condensates of the NS-NS form 
field (H 2 ) or R-R form fields of extra dimensions inte- 
grated over the whole extra dimensions 



— n ~ / dx 
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In the following we treat A as a free parameter which can 
be either positive or negative. 

The action, and thus also field equations, can be writ- 
ten in a more familiar form in the Einstein metric, which 
is defined by 



9tiu = e 2 ^~gpv (6) 

We will work in this metric from now on. In terms of 
Einstein metric, the action becomes 

1 
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The field equations derived from the action J7J are 
B p „ - \e-^B% (BB*) 



-m 2 B e 3 * 



1 + ie- 4 *£ 2 - ie" 8 * (BB*) 
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where the dilaton potential is 
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and the energy-momentum tensor is given by 

= -^Ae 2 * + 2V Al $V,$-.g^(V$) 2 



(11) 
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In the subsequent sections, we examine the equations 
of motion ll81)- (|10fl and find cosmological homogeneous 
solutions. Cosmological implications of the D-brane 
world is of our main interest, including the stabilization 
of the dilaton and the evolution of the 2-form field. 



III. COSMOLOGICAL HOMOGENEOUS 
SOLUTIONS 

A. The equations of motion 

In this section we study cosmological solutions of the 
D-brane universe in the presence of both the dilaton and 
the NS-NS form field. We assume spatially homogeneous 
configurations for the NS-NS form field and the dilaton, 
and look for the time evolution of these fields and the 
expansion of the universe. 

The non- vanishing homogeneous antisymmetric tensor 
field, in general, implies the anisotropic universe. To con- 
sider the simplest form of anisotropic cosmology, we take 
only a single magnetic component of B^„ to be nonzero, 
namely B 12 (t) = B(t) and B oi (t) = B 23 {t) = B 31 {t) = 0. 
Then the metric consistent with this choice of field con- 
figuration is of Bianchi type I 



ds 2 = -dt 2 
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Then the field equation for B, Eq. JSJ) becomes 



m 2 e^B 

- 0, 
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the dilaton-field equation © is 
'hi o 2 a 3 
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and Einstein equations IjlOjl are 
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= £>* + QB + Qb + QA, (16) 

= -£* + £>B + £b + £>A, (17) 

= -£>$ + £ B + O b + £>A, (18) 

= -£* - B + £> b + £>A, (19) 



where 
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When raj =/= 0, it is more convenient to employ the 
dimensionless time variable t — mst and to introduce 
the variables and b defined by 
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= hl( 
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B 
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Then the equations of motion are written as 
b + (di + d 2 + a 3 )b 



+ 



2$ + 2(-ai - d 2 + d 3 - 2$)$ 



+ di + d 2 + (di + d 2 )d 3 
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& = 0, (22) 



$ + (di + d 2 + d 3 )$ = -2p B - p A - TjPb ~ Z 5 ^ ( 23 ) 



aia 2 + a 2 a 3 + a 3 ai 
di + di(di + d 2 + d 3 ) 
d 2 + d 2 (di + d 2 + d 3 ) 
d 3 + d 3 (di + d 2 + d 3 ) 



p® + Pb + PA + Pb, (24) 
PA + Pb, (25) 
PA + Pb, (26) 
2pB+pA + Pb, (27) 
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where 



PA = Ae 
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where A = A/m%. In the subsequent subsections, we look 
for the solutions to the above field equations for various 
cases beginning with some simple solutions. 



B. The massless limit (tub = 0) 

To see the effect of the brane on the spacetime dynam- 
ics, let us first consider the limit of mg = which corre- 
sponds to either the absence of the brane or the limit of 
vanishing string coupling <? s , namely the usual massless 
antisymmetric tensor field. In this limit, the equation 
(|14J) is easily integrated to yield a constant of motion 



a 3 B 
aia 2 



L 3 (constant). 



(29) 



With the vanishing potential, B manifests itself by non- 
vanishing time derivatives. In the dual variable, it corre- 
sponds to the homogeneous gradient along a; 3 -direction. 
The spacetime evolution with the dilaton rolling in this 
case was studied in Ref . |Tl[ • Here we have assumed that 
the dilaton is stabilized by some mechanism. Following 
Ref. 1 1 1 1] , we introduce a new time coordinate r\ via the re- 
lation drj = L 3 dtj 'a\a 2 a 3 . Then the equations (|l(j[l 1191) 
can be written as 



*l a 2 + a 2 a 3 + a 3 a 'l = 



1 2 2 
-0^2, 



a-t = a 2 



= 0, 



2 2 



(30) 
(31) 
(32) 



where the prime denotes the differentiation with respect 
to r]. 

The solutions for a± and a 2 are trivial 



Oil = Cir), a 2 = C 2 r], 



(33) 



where Ci )2 are constants and we omitted the integration 
constants corresponding simply to re-scaling of scale fac- 
tors. The (^-equation ll-'il't is also easily integrated to 



give 
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Wi + C2) 2 
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The constraint equation (fTTTTf) relates C\^ 2 and C 3 by C3 = 
—CiC 2 j{C\ -I- C2). Then the relation between r\ and L^t 
is explicitly given by 



L 3 t = 



drj 01(77)02(77)03(7?) 



dr\ exp 
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where x = e 2{Cl+C2 > r > /8(d + C 2 ) 2 and P = (d + C 2 - 
Cs)/2(Ci + C 2 ). The evolution of scale factors for large 
L 3 t is given by 



«1,2 « (log L 3 t) q 



o 3 (x L 3 rj, 



(36) 



where qi. 2 — C\^ 2 j2(C\ + C 2 ). Therefore, with non- 
vanishing B\ 2 (t) 1 only 03 grows significantly and the spa- 
tial anisotropy develops. This can be seen clearly by 
considering the ratio H 3 / B\ 2 where Hi = di/a*. Taking 
C\ = C 2 and thus gi )2 = 1/4, we get 



Hi, 



Hog L 3 t 



(37) 



which grows as time elapses. Finally, we remark that 
such anisotropy cannot be overcome by some other type 
of the isotropic energy density in the expanding universe. 
Assume that the isotropic universe (a^ = a) is driven 
by, e.g, radiation energy density pr. Then, one finds 
Pb oc 1/a 2 from 1(20) 1 and ((20) and thus pb/pr c< a 2 , 
which implies that the late-time isotropic solution can 
be realized only in a contracting universe \H\ . 



C. The B oscillation 

Let us now take into consideration the effect of space- 
filling D-brane. To get sensible solution, we fine-tune 
the bulk cosmological constant term to cancel the brane 
tension, that is A = —m 2 B /2, so that the effective four- 
dimensional cosmological constant vanishes. We assume 
again the dilaton is stabilized in some way and set a\ = 
a 2 . Then we can rewrite the full equations as follows 



a 1 + 2a\a 3 

(ii + di (2di + a 3 ) 
a 3 + a 3 (2di + a 3 ) 
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(38) 
(39) 

(40) 
0. (41) 



First we examine the evolution of b(t) and scale factors 
qualitatively. Suppose b starts to roll from an initial value 
bo, while the universe is isotropic in the sense that dio = 
d3Q. We assume initially aio = 030 = 1 (ctio = a 30 = 0) 
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and B = so that b — B and 6 — ~2ai B . While 
b is much larger than unity, the rapid expansion of ot\ 
due to the large potential proportional to m 2 B b drives b 
in feedback to drop very quickly to a small value of order 
one. A numerical analysis shows that this happens within 
mst < 2 up to reasonably large value of bo for which the 
numerical solution is working. The behavior of b(t) after 
this point is almost universal irrespective of the initial 
value bg if it is much larger than unity. 

Once b becomes smaller than unity, the quadratic term 
of mass dominates over the expansion and b begins to 
oscillate about b — 0. Then the expansion of the uni- 
verse provides the slow decrease of the oscillation am- 
plitude. The situation is the same as that of the coher- 
ently oscillating scalar field such as the axion or the mod- 
uli in the expanding universe. For small 6, the energy- 
momentum tensor of the oscillating B field is given by 
Tf/ 1 = diag[-/3,pi,p 2 ,P3] where 



of the solutions are universal irrespective of initial con- 
ditions. This is an attractive property of the exponential 
potential. For the potential of the form V(Q) = Vbe' 3 *, 
there is an attractor solution 



$ = In 





(3V 1/2 t 



a=[- \ lnt, 



(46) 



for < (3 < vl2. The scale factor obeys the power-law 
time-dependence, implying that the rolling of $ consti- 
tutes the matter having an equation of state p = wp 
where w — (3 2 /6 — 1 varies from —1 to +1 for the afore- 
mentioned range of (3. 

We found this type of particular solutions of the 
Eqs. H22 [) -H27 p . which can be found when we have a single 
exponential term in the potential, that is, for the case of 
A = and for the case of raj =0. For both cases we 
start from an ansatz of the form 



p = \(b + 2di&) + \m\ ( \A + b 2 -l 
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1 I: 
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^(b 2 + m 2 B b 2 ), 



(42) 
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Pi = P2 = — |6 + 2di& mt , 

4 V J 2 b VVTTP 

- -~(* 2 -m|& 2 ), 



(43) 



P3 



l(6 + 2d 1 6) 2 -im|(VTT^-l) 

\ (b 2 m 2 B b 2 ) . (44) 



With the expansion of the universe neglected, the equa- 
tion of motion for 6, Eq. 141|) is then approximated by 



0. 



(45) 



Since the oscillation is much faster than the expansion, 
we can use the time-averaged quantities over one period 
of oscillation for the evolution of spacetime. The equation 
(|43|l gives the relation (6 2 ) = {m 2 B b 2 ). Thus, the oscillat- 
ing B field has the property Pi,P2,J>3 ~ and behaves 
like homogeneous and isotropic matter. This justifies the 
name of B-matter. Therefore, after b begins to oscillate, 
the isotropy of the universe is recovered. 



D. Attractor solutions 

Now, we consider the generic case that both the dilaton 
and the NS-NS 2-form field evolve in time. The dilaton 
in the system of equations (|22|I -H27 () has the exponential 
potential up to the correction due to the NS-NS 2-form 
field. It is well-known that the scalar field with the expo- 
nential potential possesses the scaling solution in which 
the energy density of the scalar field mimics the back- 
ground fluid energy density 0, 0] . This scaling solu- 
tion is also an attractor, so that the late time behavior 



ai = 7i hit, a 3 = 73 lnt, 
$(*) = 7$ In* + $o, b = constant, (47) 

where we suppressed constant terms in a\ and a 3 which 
correspond to simple rescaling of coordinates. For the 
case of A = 0, we obtain two distinguished solutions 



and 



ai (t) = a 3 (t) = -lnf, 
y 

$(t) = -H]nt + l n H ) 6 = 0, 



10 - 3 - 

ai(t) = — lnt, a 3 {t) = -lnf, 



(48) 



2 - 2, /8-5 1 / 4 
-lnt + - In 

3 3 V 21 



6 = ±-. (49) 



The first solution is nothing but the solution (|46|l with 
(3 = 3. The second solution has the non- vanishing NS-NS 
form field. For the case of = 0, we introduce a new 
rescaled dimensionless time variable t = A 1 / 2 /: instead 
of t and then get the continuous set of solutions from 
Eqs. 

ai(t) = »3(f) = lnf, 
$(t) = — lnf + — In 2, b — arbitrary constant. (50) 

<&(f) and a(t) are same as those in Eq. (|46|) with (3 = 2, 
while we have the non-vanishing B field condensate. 

These solutions are the solutions to the Eqs. I|22|l - H27(l 
for the specific initial conditions. However, the impor- 
tance of these solutions, as noted in the paragraph above, 
arises from the fact that they are attractors, which means 
that after enough time the solutions with different initial 
conditions approach these solutions. We will confirm this 
through numerical analysis in the next subsection. 

The solution (|48|) applies for the brane tension domi- 
nated case where the dilaton potential is approximated 
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by V($>) = ^m^e 3 *. The evolution of the dilaton under 
this potential produces matter with the equation of state 
p=\p. 

Once the antisymmetric tensor field is turned on, the 
anisotropy appears as in the solution l|49|) . The measure 
of anisotropy is 



This result is contrasted with that in Ref. || where the 
dilaton is assumed to be stabilized. The rolling of the 
dilaton makes the difference. It affects the dynamics of 
B field in such a way that b remains constant at b = ±1/2 
instead of oscillating about the potential minimum 6 = 
and the anisotropy is maintained. 

Let us turn to the third solution JSUJ. It is relevant 
when the dilaton potential arising from the curvature of 
extra dimensions, V($>) — Ae 2 *, dominates over other 
contributions. In our scheme this happens as the dilaton 
rolls down the potential. When b vanishes, the transition 
point at which Ae 2 * starts to dominate over im^e 3 * 
is at $t = ln(2A). Thus, this solution describes the late 
time behavior of all the solutions with various initial con- 
ditions when A is positive. It is very intriguing since we 
achieve the isotropic universe in the end. The anisotropy 
is the result of the non- vanishing B field, whose coupling 
is dictated by the brane tension term. The rolling of 
dilaton makes the brane tension term less important than 
the extra dimensional curvature term, which recovers the 
isotropy. The rolling of the dilaton under the potential 
Ae 2 * now forms the matter having the equation of state 
p = — ip, thus giving marginal inflation. 

E. Numerical analysis 

1. Initial conditions 

We have the second order differential equations for four 
variables B(t), a\{t) — (^(i), a^{t). Thus we need 

eight initial values $ , $0, B , B , a w , "io, «30, a 30 
to specify the solution. Among these, c^o can always be 
set to zero by coordinate rescaling. c^o must obey the 
constraint equation l|24ll . but this does not fix the ratio 
dio/a3o. We choose the isotropic universe as a natural 
initial condition which leads to dio = 0:30- 

Since the dilaton potential is composed of exponential 
terms, the shift of the dilaton field by a constant can 
be traded for the redefinition of mass scale. We use this 
property to take the initial value of the dilaton to be zero 
without loss of generality. In our numerical analysis, we 
take the dimensionless time variable as t = rht where 
m = TOse^* and use the variable $ = $ — $0 with 
its initial value $0 — 0. This means that the proper 
time scale for the cosmological evolution is not m^ 1 , but 
mT 1 . For convenience's sake, we take the initial time as 
to = 1. The other mass scale A is also affected by this 



shift and we can treat it by replacing the parameter A 
withA=Ae-*«. 

Now we need three initial values $0, Bo, and Bo, to 
fix the functional form of the solution. The initial values 
60 and 60 are related to Bo and Bq by 60 = Bo and 
b = B - 2(di + $o)5o- 

2. Numerical solutions for B 7^ 

We turn on the NS-NS 2-form field along the a; 3 direc- 
tion, #12 = B =/= 0. The spacetime becomes anisotropic, 
ai(t) — ci2(t) ^ 03(4), in general. The solutions are clas- 
sified by the signature of A. For A < 0, the solution 
becomes singular. Here we skip the description of such 
singular solutions which are not suitable for the evolution 
of our Universe. 

For A > 0, the evolution is divided into two stages. 
In the first stage where pb is dominant, the solution ap- 
proaches an attractor l|49l) of the A = case. The ap- 
proached value of b is either +~ or — ~ depending on the 
initial conditions. If we look at the evolution of each 
component of energy density, the kinetic energy of the 
dilaton p$ catches up the potential energy pb and the 
ratio of them becomes constant. This is a characteristic 
feature of the scaling solution ^|. The kinetic energy 
of B field is kept much smaller than both of them, but 
the anisotropy is still maintained due to the difference 
between pb and pb- In the second stage where pA is dom- 
inant, it approaches another attractor JSUJ. Thus the 
universe recovers the isotropy. The final value of b is a 
certain constant which is determined by initial conditions 
and can differ from ±i. Numerical solutions for a few 
initial conditions are shown in Figure ^ The kinetic en- 
ergy of the dilaton /?$ catches up the potential energy pb 
in the first stage and pa in the second stage. The ratios, 
p<t,/pb and p$/pa, approach constants in each stage. The 
kinetic energy of B field is kept much smaller as in A = 
case. 



IV. THE DILATON STABILIZATION 

The vacuum expectation value of the dilaton deter- 
mines both the gauge and gravitational coupling con- 
stants of the low energy effective theory. Therefore, the 
dilaton must be stabilized at some stage of the evolution 
for the action Q to have something to do with the real- 
ity. In this section, we study the cosmological evolution 
when the dilaton is stabilized. As for the correct mecha- 
nism of dilaton stabilization, the consensus has not been 
made yet. Our goal here is to illustrate an example of 
the dilaton stabilization and look into the effect of it on 
the dynamics of the NS-NS 2-form field and the cosmo- 
logical evolution, since the overall features of which are 
insensitive to the detailed mechanism of stabilization. 

Our starting point is the dilaton potential (|11|) . This 
potential possesses the minimum for A < 0, but the value 
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FIG. 1: Numerical solutions for B 7^ and A > 0. When A = 10 4 , the solutions of four initial conditions are given: 
$ = -7/4 ■ 5 1/4 , 6 = 1/2, B = for the solid curves, $ = 0, b = 1, B = for the dotted curves, $ = 0, b = 10, B = 
for the dashed curves, and $0 = — \/3/2, 60 = 1, Bo = for the dash-dotted curves. The lower-right panel shows the evolution 
of each component of energy density for initial conditions $0 = — \/3/2, 60 = 1, -Bo = 0. 



of the potential at the minimum is negative and needs to 
be set to zero by fine-tuning of the constant shift. How- 
ever, the constant shift of the potential has no motivation 
in the context of string theory. Instead, we introduce a 
term ^m|e~ 3 * in the potential, which can arise from the 
effect of various form field fluxes in extra dimensions ^3] • 
Thus, our potential for the dilaton for B = looks like 

^ F (d>) = I ( TO | e 3* + 2Ae 2 * + m 2 F e" 3 *) . (52) 

This potential has a global minimum for any value of A 
and mp. To have sensible cosmology, the potential at 
the minimum must be zero. For A < 0, this can be done 



through a fine-tuning of the parameters in the potential 

" 2 =K-s A ) 6 ' < 53 > 

where /1 2 = m\ /m 2 B . Then the minimum is located at 
$ F =m(-^ . (54) 

The shape of this fine-tuned potential for A = —0.1 is 
shown in Figure |2J (a). Now the equations l|23|) - (|27() are 
modified accordingly 
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$ + (cVi + a2 + a 3 )<I> = -2p B - pa- -Pb- Pb + 7;PF, (55) 




(c) (d) 




i t 

FIG. 2: (a) The potential V{<&) for A = —0.1. The minimum of the potential is at $ m = — In 6 with V($ m ) = 0. (b), (c) and 
(d) show the numerical solutions of <&(t), 6(1) and au(t), respectively, for A = —0.1 and initial conditions <3>o = 0, $o = and 
bo = 1. $(i), b(t), and tdti(t) approach $ m , 0, and 2/3, respectively. 



&i + ai(di + d 2 + A3) = PA + Pb + PF, (56) 
d 3 + d 3 (di + d 2 + d 3 ) = 2p B + PA + Pb + pF, (57) 

where p-p = ^p 2 e~ 3 *, while the equation (122[l for 6 is not 
changed. 

With the stabilizing potential for the dilaton, the cos- 
mological evolution is completely changed. The dilaton 
and the NS-NS 2-form field rapidly come to the oscilla- 
tion about the potential minimum $ = <5>p and 6 = 0. 
Oscillating $ and B fields behave like ordinary matter 
satisfying the equation of state p = ||. The universe 
becomes isotropic and matter dominated. The numer- 
ical solutions for the stabilizing dilaton potential l|52|) 
with A = —0.1 are plotted in Figure IU-(b) to (d). One 



can see the oscillation of $ and b, and that both di and 
d 3 approach to 2/3t, indicating the matter domination. 
Damping of $ and b oscillations is due to the expansion 
of the universe. 



V. CONCLUSION 

We investigated cosmology of a four-dimensional low 
energy effective theory arising from the NS-NS sector of 
string theory with a D-brane which contains the dynam- 
ical degrees of freedom such as the gravity, the dilaton, 
and the antisymmetric tensor field of second rank, cou- 
pling to the gauge field strength living on the brane. The 
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dilaton gains a potential in the presence of the D-brane, 
the fluxes of various form fields and the curvature of extra 
dimensions. The NS-NS 2-form field becomes massive in 
the presence of the D-brane. The dynamics of the system 
crucially depends on the curvature parameter A and the 
brane tension parameter ms through which the dilaton 
obtains the potential of the form; Ae 2 * + im^e 3 *. Here, 
the latter becomes the effective mass of the 2-form field. 

When the 12-component B(t) of the 2-form field is 
turned on, the universe undergoes an anisotropic ex- 
pansion described by the Bianchi type-I cosmology. We 
found the attractor solutions showing the overall features 
of general solutions and confirmed it through numerical 
analysis. The dilaton $(t) runs to the negative infin- 
ity settling to a logarithmic decrease in time. When the 
brane tension term ^m^e 3 * dominates, the anisotropy 
due to the 2-form field flux is sustained. If there is a pos- 
itive curvature term Ae 2 *, it dominates finally over the 
brane tension term as the dilaton rolls down to the neg- 
ative infinity. Then the expansion of the universe turns 
to be isotropic and linear in time. Accordingly, B(t) de- 
creases inversely proportional to time. 

For sensible phenomenology and cosmology, the dila- 
ton must be stabilized. In order to study the dynamics 
of the 2-form field (B-matter) and the stabilized dila- 
ton system, we adopted a dilaton mass term of the form 
m|e -3 *. Then the dilaton potential has a global mini- 
mum and the cosmological constant of our Universe can 
be fine-tuned to a desired value with negative A. With 



this stabilizing potential, we obtain a reasonable cosmol- 
ogy from an initially anisotropic universe resulting from 
the 2-form field flux. While the dilaton evolved to a sta- 
bilized value, B(t) shows an oscillatory matter-like be- 
havior (B-matter), and the universe expands as in the 
usual matter-dominated era recovering the isotropy. 

Finally, there are other dynamical degrees of freedom 
in the D-brane world, such as the R-R form fields, which 
are not included in the present work. Investigating the 
cosmological implications of them is left for the future 
work. 
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